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1. DESCRIPTION OF STUDIES UNDERTAKEN 


1 1 Scope 

The objective of the studies in this report was to develop rotordynamic 
analysis procedures for obtaining numerical criteria for the balancing 
of flexible rotors. Present numerical criteria for flexible rotor 
balancing are based on related practical experience. This study 
attempts to develop flexible rotor balancing criteria, through the use 
of rotor analysis techniques, Such a formulation should validate exist- 
ing criteria, and would provide a general procedure for establishing 
criteria values for specific rotor types. 

The procedure used in this study was to use a rotor response computer 
program to calculate the response of a range of flexible rotors in a 
systematic manner. This provided rotor response charts by which criteria 
for many similar rotor types could be obtained directly, as follows. 

First, a comprehensive analysis was made of the unbalance response of 
a uniform flexible rotor in several types of fluid-film bearings. 

Charts 6f rotor whirl amplitude vs. speed were plotted for ranges of bearing- 
rotor stiffness parameter, C/6. Second, a procedure for representing a 
non-uniform rotor as an equivalent uniform rotor with similar dynamic 
properties was developed. This was based on using the lowest rigid- 
bearing flexural critical speed value and the mass of the rotor to de- 
fine an equivalent rotor stiffness, which was then matched to that of 
the uniform diameter rotor. A further procedure was developed for mak- 
ing equivalent bearing representations. This is important in cases where 
a rotor is supported in two end bearings which do not have identical stiff- 
ness and damping properties. 

The third step was to test the equivalent rotor procedure for its 
ability to represent bearing whirl amplitude,;. This was done by 
computing the response of the equivalent rotor system and that 
of the actual non-uniform rotor system, and comparing results. This 
comparison was made for twelve different non-uniform rotor cases, to 
compare the effectiveness of the equivalent rotor method amplitudes for 
the first whirl mode. 






The fourth step in the procedure was to compare results from the 
equivalent rotor method with results available from other calculations 
done on six practical rotors. In several cases test data on rotor 
whirl amplitudes was available for these practical rotors. This data 
allowed information on rotor unbalance/rotor journal whirl amplitude 
and on other response details to be developed. Equivalent .’otors were 
developed in each instance, and whirl amplitudes were calculated for 
specified rotor residual unbalance. These results were checked against 
results based on actual (-^) values. The results provided an assess- 
ment of the validity of the proposed method, and gave new information 
on the suitability of existing rotor balance criteria. 

The above four steps are demonstrated herein, and typical results are 
shown. The remainder of this report discusses progress made toward 
improved flexible rotor criteria, with suggestions for additional studies 
to supplement the above. 

1 . 2 Need 

The need to develop improved balancing criteria arises from both 
the lack of an established data base for flexible rotor balancing, 
and from the lack of a procedure by which specific criteria for a 
given rotor can be obtained. Such data have been developed for rigid 
rotors in two end bearings. The method for selection of an appropriate 
numerical criterion for a rigid industrial rotor requires that the rotor 
be classifie into a specific quality grade suited to its function. 

A corresponding residual unbalance value is then selected for its 
speed of operation. Charts for both the above steps are given in 
I.S.O. Standard Document 1940 (1973): Balancing of Rigid Rotors, 

reference [1]. The charts in this document are based on comprehensive 
data for rigid rotors from a range of industrial applications: see 

Muster and Flores [2]. 

Similar data for flexible rotors are in the process of practical de- 
velopment. This development is presently based on two sources: (a) 

from established machinery vibration standards and, (b) from prior 




experience with typical flexible rotor machines. Data obtained from 
both procedures is empirical, Without questioning its validity, the 
data base for such information must be somewhat restricted. Also, it is 
difficult to quantify in a consistent manner. These reasons suggest 
that additional data developed by analysis should represent a useful 
addition to existing criteria, Once developed, an analytical criterion 
would apply to a rotor of any shape, size, or speed. 

1.3 Previous Work 

Early computer studies of unbalance response of a rotor in fluid-film 
bearings were made by Lund and Sternlicht [3] in 1961. A single-disk 
rotor was studied to deteniiine the influence of the oil film in attenua- 
ting rotor vibrations, especially during operation near the flexural 
critical speed of the system. Data on the dynamic properties of several 
types of fluid-film bearings is included in this report, with comparisons 
of their relative influence on rotor unbalance amplitudes. The response 
of a symmetrical two-disk rotor in partial bearings was studied by Warner 
and Thoman [4], who gave design charts for evaluation of amplitude re- 
sponse and transmitted force. The computer analysis of rotor-bearing 
systems has been discussed by several authors including Lund and Orcutt 
[5] who compared theoretical whirl amplitudes with test results for a 
given high-speed flexible rotor. These results were later verified by 
Thomas [6] who wrote the initial version of the computer program described 
herein. A parametric study of the unbalance response of a uniform flex- 
ible rotor in fluid-film bearings was made by Rieger [7]. 

Criteria for rotor balancing were initially proposed by Feldman [8] 
based on work by Federn [9], Rathbone [10], and others. Muster and Flores 
[2] collected data on rigid rotors from a variety of industries, and from 
this data made a statistical determination of acceptable rotor vibration 
levels for various rotor types, sizes, weights and speeds. These results 




were used as the basis for the rigid rotor criteria charts in I.S.O, 
Document 1940 (1973) [1]. This document also contains much valuable 
guidance on rigid rotor balancing technology, Other standards for 
rigid rotor balancing (A6MA, NEMA, API) are also available, but I.S.O. 
1940 is more comprehensive. 

Flexible rotor balancing procedures and concepts are described in 
I.S.O. Draft International Standard I.S.O. DIS 5343 (1976), ref. [11], 
This document describe;, suitable procedures for balancing the five basic 
rotor classes mentioned previously. Numercial criteria for flexible 
rotor balancing were proposed by Giers [12] based on values given in a 
German Vibration Standard Document VDI 2953. In revised form these 
criteria plus a guide to their application is given in I.S.O. draft 
International Standard 2905 (1978), ref. [13]. These criteria values 
make allowance for rotor size, quality, operating speed, and location 
of measurement. Further detail form I.S.O. 5343 and I.S.O. 2916 are 
given in the present report. 




2. UNBALANCE RESPONSE COMPUTER PROGRAM 


2.1 Program Details 

The computer program used in this study was a general rotor system 
unbalance response program, based on the system model shown in figure 1. 
Practical rotors are modeled using a sequence of cylindrical sections 
which are rigidly joined together at their ends, Disk sections of the 
rotor with mass and gyroscopic properties are also incorporated at 
these end junctions, as shown. Rotor unbalance is included at the 
rotor junctions, and a spatial distribution of unbalance along the rotor 
may be introduced through the use of sine and cosine unbalance components 
at the rotor unbalance stations. 

Bearings are represented in this program by the model shown in figure 2, 
using the conventional linearized eight stiffness and damping coefficient 
procedure. Bearings are introduced into the system model at suitable 
rotor section ends, as shown in figure 1. This procedure allows the 
influence of any desired bearing geometry or bearing operating conditions 
to be investigated for which the coefficient data is available. The 
influence of given plain cylindrical bearings, tilting-pad bearings, 
axial groove bearings, and partial arc bearings on the response of 
several rotor configurations was investigated during this investigation. 

The rotor-bearing formulation used allows the shaft to whirl in a elliptical 
orbit, where required by the bearing conditions. The magnitude of the 
rotor whirl ellipse maximum radius is plotted for selected rotor locations 
on the rotor response charts. These response charts are used in developing 
the balance criteria with which this report is concerned. 

The computer program was developed and verified by Thomas [6]. A matrix 
formulation is used, based on the dynamic stiffness procedure of 
McCallion and Rieger [14]. Excellent correlation with previous rotor- 
bearing system studies was demonstrated during program verification. 

A program listing with input and output details is given in the Appendix 
of this report. 




3. UNIFORM ROTOR RESPONSE CALCULATIONS 


3.1 System Parameters 

It has been shown by Rieger [7] that the unbalance response of a uniform 
rotor in fluid-film bearings may be described in terms of: 

a) bearing geometry and operating conditions 

b) axial location of unbalance 

c) type of unbalance (force, couple) 

The dimensionless parameters used to describe rotor unbalance response 
are: 


1. Amplitude ratio 


II - maximum radius of whirl orbit, ins, U = Wa 

a " eccentricity of rotor mass, ins. 


2. Speed ratio 

w _ shaft speed of rotation, rad/sec 

to ” lowest rigid bearing critical speed, rad/scc 

3. Stiffness ratio 


c_ _ radial machined clearance of bearing 

6 “ Mid-span deflection of uniform rotor in rigid bearings 


4. Eccentricity ratio 


£ _ eccentricity of journal center from bearing center 
c ~ radial machined clearance of bearing 




3 • 2 Uniform Rotor Response Charts 

Figures 4 through 29 show charts of unbalance response for a uniform rotor 
In: 

a) Plain cylindrical journal bearings 

b) Four-shoe tilting pad bearings 

c) Partial arc bearings 

d) Four axial groove bearings 

Details of the conditions for which the above results apply are specified 
in table 1, These results apply fcr a uniform cross-section end-bearing 
rotor as shown In figure 3. In each instance, the unbalance force is a 
single force vector located at the axial position specified in table 1. 

The unbalance response computer program described in section 2 was used 
to calculate the rotor response results. In most Instances, the rotor 
was represented in two sections as shown in figure 3. These sections 
join at the unbalance vector location. This representation includes the 
distributed mass and distributed stiffness properties of the rotor, which 
are Incorporated directly through the rotor dynamic stiffness mentioned 
previously. The other ends of the rotor sections are each supported In 
fluid-film bearings of the type specified. Ti.e validity of results 
obtained using this rotor model was verified against data from ref [7] 
before proceeding. 

Data given in the charts of unbalance response applied for conditions of 
constant bearing eccentricity ratio, e. This representation simplified 
the calculation procedure. To use these charts it is necessary to know 
in advance the bearing operating eccentricity for the conditions under 
which the rotor unbalance response data is required. Bearing eccentricity 
can be obtained from standard steady-state performance charts, for the 
appropriate bearing type, slenderness ratio (L/D), and Sommerfeld number 
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TABLE 1. Details of Uniform Rotor Cases Calculated (Con*t} 
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for the given bearing operating conditions, Bearing eccentricity may also 
be obtained by a measurement from an actual bearing. Additional details 
of pad length, groove geometry, etc., for which the bearing coefficients 
apply must also be considered. When applying the results of this report 
it should be borne in mind that variations from the calculated bearing 
geometry details and operating conditions may cause response results to 
differ somewhat from the computed results, especially In the region of 
the critical speed peak amplitudes. 

3.3 Uniform Rotor Response Results 

Figures 4 through 29 are charts of rotor whirl amplitude results 
obtained using the computer program described in Section 2. These 
results indicate the variation of: 

a) bearing type (plain, axial groove, partial arc, tilti' pad) 

b) stiffness ratio (c/6 « 0.3, 1.0, 3.0, 10.0, 30.0) 

c) speed ratio (w/wn = 0,1 through 25.0, relative to rigid bearing 
critical speed, Wf,) 

d) response location (midspan response, journal response) 

The curves shown apply for the case of midspan unbalance. The charts 
are discussed in detail below. 

a) Plain cy l indrical bearings 

Figure 4 shows the variation of mid-span rotor response with. mid-span 
unbalance over the speed range at eccentricity ratio c = 0.2, for various 
stiffness ratio (c/6) values. Low stiffness ratios, e.g., c/6 =0.3 
correspond to the rigid bearing case (high 6, low C), which is seen to 
peak at w/wn = 1-0, and 9.0, as expected. The peak at w/wn = 4.0 is 
absent because the mid-span unbalance does not excite non-symmetric 
modes. Higher stiffness ratio values tend toward rigid-rotor flexible- 
bearing conditions. For c/6 = 0,3 the minor peak at w/w^ = 0.5 is a 




clamped rigi'd-roter critical speed, and the peaks at 2.25 and 12.0 
represent free-free rotor modes. At u/wn = 2.25 the bearing damping 
is still significant, whereas at 12.0 the bearing damping is less effective 
and large whirl amplitudes occur. Similar trends are evident in figure 5 
for B 0.5. More rigid bearing conditions exist with e = 0.7 as shown 
in figure 6. 

For this case, the rigid bearing curve c/d = 30 shows two low-speed 
peaks at - 0.2 and 0.5, corresponding to the damped rigid rotor 
critical speeds. This stiff bearing case also brings out the free-free 
criticals at w/w^ = 2.25 and 12.0 (strongly). 

Results for journal whirl amplitudes with mid-span unbalance are shown 
in figures 7, 8 and 9, corresponding to eccentricity ratios of e = 0.2, 

0.5 and 0.7. Amplitudes are generally less than at mid-span, though 
similar trends may be observed in the results. For e - 0,2 the rigid 
bearing case c/6 0.3 shows resonant peaks at w/wp = 1.0 and at 9.0, as 
previously. The rigid rotor case shows resonant peaks at w/wp = 2.25 
and 12.0 as previously, but the low-speed rigid rotor modes below 
w/a)f| = 1.0 are evidently suppressed probably by over-critical damping 
in the bearings. Similar results are apparent for c = 0.5. The stiff 
bearing case e = 0.7 shows the rigid body critical speeds below = 1.0, 
possibly because the bearing dam, ping is effectively smaller for this 
case. 

b) Four axial groove bearings 

Similar trends are evident in the axial groove bearing results to those 
observed for the cylindrical bearing rotor results. For e = 0.2 in 
figure 10 the rigid bearing curve c/6 =0.3 again shows a sharp resonant 
peak at w/wp = 9.0~10.0, though this appears to have been suppressed by 
damping, as was the peak around co/tO|,| = 16.0. Thus even at mid-span the 
rigid bearing rotor tends to be adequately damped in its higher modes. 

The rigid rotor trends differ from these results. A rigid body mode 
is observed at w/wp = 0.46, and the rotor free-free mode occurs at 2.25. 

The pe.ik amplitude in this mode is larger than the rigid-body mode, and 


the next res‘»nant mode at w/u)p, « U’.O is even stronger, showing the 
decreasing influence of damping for higher rotor frequency ratios. 

As the rotor stiffness ratio c/iS increases the above trends become 
more pronounced, 

For r - 0.5 figure 11 all resonant peaks are very distinct, i.e., the 
system is losing its effective damping. The rigid bearing case c/d = 0,3 
is evident at w/w,i 1.0, 9.0 and 16.0. The rigid rotor case is also 
sharp and distinct at = 0.45, 2.25, and 12.0 as discussed previously. 
All resonant peaks apart from c/d = 0.3 blend into a single response 
above m/Wj^ ”9.0 with a strong common response peak at 12.0 

For r 0.7, figure 12, the rigid bearing rotor resonances occur at 
w/w,-, = 1.0 (strong), 4.5 (damped, weak), 7,0, and 12.0 (very strong). 

The rigid rotor criticals for c/d = 30 occur at u)/w,-| = 1.5, 0.4 (damped), 
2,25 (strong) and 12.0 (very strong), c/d cu'-'ves except 0,3 appear to 
blend around a,wn = 2.0 and so form common resonant peaks at 2.25 and 
12 . 0 . 


The journal response at the above conditions is shown in figures 13 
through 15. Figure 13 shows the journal response for c =0,2 which is 
similar to that observed for the plain cylindrical bearings. The rigid 
bearing case has the lowest journal whirl amplitudes with resonant 
conditions at oi/wp = 1.0, 7,5 and 16.0. Amplitudes generally increase 
as rotor stiffness increases for c/d = 30. A rigid critical is observed 
at 3.0 (stronger) and 12.0. Again the stiffen rotor responses appear 
to form a conwiion peak around w/wp = 12.0, though this trend is apparent 
in the lower resonance at w/wp = 2,25 to 4.0. The peaks become 
sharper as the rotor stiffness increases. 

For an eccentricity of r =0.5 the resonant peaks become sharper, due 
to the overall decrease in bearing damping, i.e., (Bw/k) (effective) 
decreases. This effect is evident in figure 12, The tendency toward 
blending for stiffness ratios above c/d = 0.3 is again observed at 
w/w,^ = 8.0 and beyond. It is also evident that all response peaks have 
maximum values of about the same value, in each resonant range. 




An eccentricity of e - 0,7 the first mocle peaks are relatively small 
and of similar maximum value, X/a 6.0, as shown In figure 15. The 
second resonance at w/Wf, = 2.25 shows all stiffness characteristics 
blended, and again at w/Wf, - 12.0, with the exception of c/fi - 0.3. 
Generally speaking, the axial groove oeaks are greater than the plain 
cylindrical maxima by between 2:1 and 3:1. 

c ) Partial Arc Bearings 

Trends observed for the two previous bearing types are again evident 
for partial arc bearings. For a flexible shaft in rigid bearings and 
an eccentricity ratio of c « 0.2, figure 16, mid-span amplitude peaks 
are evident at w/wj^ - 1,0, 9.0 (minor) and oj/w^ - 12.0. Rigid rotor 
conditions give rise to larger amplitudes below oi/wn = 1.0 with damped 
rigid-body criticals in this region. Free-free modes occur at 

= 2.25 and 13.0, Again the rigid shaft characteristics begin to 
blend above w/wn = 2.0, with a common peak at w/wp = 12.0. In general 
mid-span whirl amplitudes are somewhat higher than with plain cylindrical 
bearings, 

At an eccentricity ratio of e = 0.5, figure 17, rigid-body critical 
speeds are evident for each shaft stiffness ratio above c/6 = 0.3, at 
speeds below co/Wf, = 1.0. Above w/wp = 1.0 the c/6 = 0.3 curve response 
again has critical speeds at w/w^ = 3.5, 9,0, and 13,0, whereas the 
rigid shaft criticals are at w/Wj^ = 2.25 and 13.0. Similar trends are 
evident in the curves for e = 0.7, figure 18, except that the rigid 
body criticals are more evident at the higher eccentricity ratio than 
for e = 0.2. At e = 0,7, the curves ai^e again similar to the responses 
at lower eccentricity ratios. The response amplitudes are somewhat 
higher than at e = 0.2 and 0.5. 

The journal whirl characteristics shown in figures 19, 20 and 21 are 
similar to those for the plain cylindrical bearing. A minor peak occurs 
for the flexible shaft case at w/wp = 1.0, for stiffness of c/6 = 0,3 
all other stiffness ratios show common resonant peaks at to/o)p and 2.25 
and 12.0. Magnitude of the peaks is similar to the plain cylindrical 
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journal amplitude ratios. At e = 0.5 the curves are all blended above 
w/w^ = 1.0, with peaks at 2.25 and 12.0. Below = 1.0, minor rigid 
rotor critical speed peaks are evident. These peaks are again evident 
for e = 0.7, and the rigid rotor peaks at 2.25 and 12.0 are again present. 
The flexible rotor characteristic shows a small rigid bearing critical at 
co/wn = 0.8 and at 3.5. The lower critical speeds suggest that the 
partial bearing acts somewhat more flexible than the plain cylindrical 
bearing for given stiffness ratio. 

d) Tilting pad bearings 

Response characteristics with tilting pad bearings shown in figures 22 
through 27 again resemble the characteristics seen with plain cylindrical 
bearings. The rigid bearing flexible shaft curve shows critical speeds 
at w/ti)p = 1.0 and 12.0. The higher critical is more heavily damped, and 
both peaks are lower than for corresponding plain cylindrical bearing 
case, in the ratio 0. 4:1.0 approximately. The indication is that tilting 
pad bearings dampen flexible shaft modes more than plain bearings. 

Rigid rotor criticals of low magnitude appear below tu/ojn = 1.0, and at 
= 2.25 and 12.0. The magnitude of these peaks is of the same order 
as with the plain cylindrical bearing. 

At e = 0.5, figure 23, the above pattern is repeated. The flexible shaft 
c/6 = 0.3 is lower at oi/tOf, = 1.0 and higher at cu/tjop = 12.0. The rigid 
rotor peaks are higher at low speed, but otherwise similar. At e = 0.7, 
figure 24, conditions are similar to c = 0.5 and the peaks are of the 
same general magnitude. 

The journal response curve for e - 0.2, figure 25, also resembles the 
response seen in plain cylindrical bearings. Flexible shafts in rigid 
bearings have lower response amplitudes below o)/wn = 1.0, and then grow 
to a peak at = 12.0. Rigid rotor response is larger, but no 
critical peak appears below = 1.0 for e = 0.2. The free- free peak 
appears at co/oJj^ = 2.25, and all curves blend closely into the peak at 
= 12.0. Similar trends are evident at e = 0.5, figure 26, and 
again at e = 0.7, figure 27. Amplitudes appear to be larger at higher 
operating eccentricities by between 25 to 35 percent over amplitudes 
at e = 0,2. 










it, THE EQUIVALENT ROTOR BEARING SYSTEH 

^ T Equivalent Rotor Calculation Procedure 


In the present context, an equivalent rotor**bearing system is a 
uniform rotor in fluid-film bearings which has similar dynamical pro- 
perties to those of the complex rotor-bearing system on which it is 
based. The equivalent system is a convenient means of representing 
the actual rotor system by which balance criteria suited to an actual 
rotor system may be obtained from the uniform rotor response charts dis 
cussed in Section 3. Both an equivalent rotor and equivalent bearings 
may be required to model a given system. The concept of an equivalent 
rotor system is shown in figure 30, The procedure for defining the 
equivalent rotor is as follows: 

Step 1: Calculate the first bending critical speed of the actual 
rotor, mounted in rigid radial supports which are located 
at the axial midpoint of the rotor journals. 


Step 2: Calculate the diameter of the equivalent uniform shaft in 
rigid bearings which has the same bending critical speed 
N ^1 as the actual rotor in rigid bearings, given the same 
material properties E and w, and the same bearing span L. 
To do this let: 



Actual 



Equivalent 


or 


9.55 


9.55 

2:39 



Elq 





1s the diameter of the equivalent rotor in rigid bearings 
which has the same first bending critical speed as the equi- 
valent rotor in rigid bearings. 

Step 3: Determine the operating eccentricity of the bearings of dia- 
meter Dg for the equivalent uniform rotor system. The theory 
for this step is given in the following section. 

The above steps define an equivalent rotor-bearing system for which the 
unbalance amplitude response properties may be determined from the charts 
given in Section 2. 

4 . 2 Theory of Equivalent Bearing Procedure 

The theory for determining the equivalent bearing eccentricity from the 
actual bearings is based on several assumptions: 

(a) The ratio of bearing radius to radial clearance (R/C) remains the 
same for the uniform rotor as for the actual rotor. 

(b) The actual bearings are sufficiently close in size and loading 
that the operating eccentricity for both bearings is approxi- 
mately the same. 

This allows the equivalent bearings to be based on the average operation 
eccentricity, or 


e 


a 





L = left 

R = right (1 ) 


For dynamically similar bearings, the assumption is made that 



U = uniform 
S = stepped 




where R is the bearing radius, C is the bearing clearance, U refers 
to the uniform rotor, and 'a* is the average value for the stepped rotor. 
Then 



as R 


U 


( 3 ) 


The maximum static deflection due to gravity of a uniform-section beam 
in rigid supports is: 


qU „ A wAL^ w 16 L^ 

® " 384 Cr" “ 384 E“^ 


Thus (C/d)*^ may now be easily calculated. 

Now, let Cg be the average bearing eccentricity of the stepped rotor, i.e. 

_ ^‘r 

‘a " ■ 2 


then since the average Sommerfeld number of the stepped rotor = f(l/e^), 
it follows that 



If we again preserve R/C, and also L/D, y, and N (where y is the viscosity, 
and N the speed in rev/sec) then 



W 

e a 

a 





where is the weight of the equivalent uniform rotor, W^, is the 
weight of the stepped rotor, is the ixe.racje diameter of the stepped 
rotor, and D is the diameter of the equivalent uniform rotor. 

Knowing (C/d)^j and Cy, the unbalance response of the stepped rotor 
can be determined by referring to the appropriate uniform rotor curve in 
this report, or it can be calculated using R0T0R2, or another unbalance 
response program, 


4 • 3 Sample Calculation of Equivalent Rotor-Bearing System 

Consider a non-uniform steel rotor of the type shewn in table 2, case 1, 
which has a bearing span of L = 100 in. and weight W = 6224.2 lb. The 
lowest rigid-support critical speed for this system is = 6400 rpm. 
Calculate the diameter of the equivalent uniform shaft having similar 
dynamic properties and the eccentricity of the equivalent plain cylindri- 
cal bearing with L/D ratio == 1..0. 


Step 1: Equivalent diameter 0^ 


Given = 6400 rpm 

Nr , 2 

and D„ y /gg 



6400 
" 2.39 


^100^2 / 0,283 

^ yT30.10®)(386.4) 


= 13.13 ins. 

C 


Step 2: Equivalent bearing eccentricity 

’e ' -ST 03-13)''= 1669.9 in'* 















11.5 


A, 


= 1.67 


= 1.63 


= 1.71 


2.01 





” fl’e ■ 3 HI.M in^ 


Rotor mid-span deflection 8 . 


5 wAL' 
’384’ ’ Tf 


5 


384 


.66)(10^) 
(30.10^)0569.9) 


1.0896 X 10 "'^ ins. 


Effective bearing machined clearance: 


A 


(^)IO X 10'^ = 


8.953 X lO"'’ ins. 


Stiffness rat'Jo: 

S = 8.953 X 10 "^ ^ 
'^G 1,0896 X 10 "^ 


If the rotor in case 1 operates with eccentricity e, = 0,2, and the 

u 

uniform rotor weight If is 

6 


= WAL = 0 . 283 ( 141 . 6 t)('iSO) = 4008.98 lb. 

6 


the equivalent bearing operating eccentricity is: 




♦ ‘ 


D n 

(■snif 


,4008. 98 w 15.0 


)^( 0 . 2 ) 


0.16 


The parameters required for determining the unbalance response of the 
equivalent rotor-bearing system are therefore: 

Stiffness ratio: (*|) - 8,217 

Eccentricity ratio: “ 0.16 

for plain cylindrical bearings of L/D 1.0 mounted 1n rigid supports. 
These values may now be used to determine the response of the given rotor 
system using the equivalent uniform rotor system response charts. 

The following section describes a range of tests which were applied to 
the above theory, to determine whether the equivalent rotor procedure 
could be used to predict the unbalance response of non-uniform rotors, 
using the uniform rotor response charts given in Section 3. The 
general result of the following study is that for the range of two-end- 
bearing rotors shown in table 2, the equivalent rotor theory predict 
resonant response amplitudes which are approximately 1.6 to 2.5 times 
as great as the true rotor response. Accepting this excess amplitude 
as a safety factor will therefore allow the equivalent rotor theory to 
be used as a basis for selecting rotor balancing criteria, This 
concept is further tested in Section 6 on several specific rotors. 




5. STEPPED ROTOR CALCULATIONS 


5,1 Roto r System Details 

The non-uniform rotor configurations shown in table 2 were selected to 
study the effectiveness of the proposed equivalent rotor method for 
unbalance response calculation, Dimensions of each rotor are given in 
table 2, Fluid-film bearings with L/D ° 1,0 were used in each case. 

An unbalance response calculation was made with each rotor in each of the 
four bearing types mentioned; Plain cylindrical bearings, partial bearings, 
tilting-pad bearings, and axial-groove bearings. The results obtained 
are shown in figures 31 through 45, 

In each case, the rotor was assumed to be a continuous flexible steel 
structure, without any deformable joints or hinges, Each rotor was sup- 
ported in two bearings located near its ends, as shown in taole 2. Effects 
arising from shaft overhangs were examined in several cases. Each rotor 
model was composed of two or three basic shaft sections, each of which had 
differing diameters and lengths. Where the bearings have different sizes, 
they have different eccentricity ratios, denoted by and respectively. 

Three calculations were made for each rotor as follows; 

(a) Response of stepped rotor to unbalance when operating in rigid 
bearings, 

(b) Response of stepped rotor in damped, flexible bearings, 

(c) Response of equivalent rotor in equivalent damped, flexible 
bearings. 

Table 2 gives lasults for the ratio of the mid-span amplitude of the equi- 
valent rotor to the mid-span amplitude of the stepped rotor, both 

in flexible bearings. It is shown that the equivalent rotor whirl ampli- 
tude exceeds the stepped rotor amplitude for all end-bearing rotors by a 






* • 


factor of between 1.63 and 2.5. The overhung rotors tn eases 5 and 7 
appear to be strongly Influenced by whether a mid-span mass exists, and by 
the overhangs theniselv*is. Further discussion of these results Is given 
In section 5.4. 

5 . 2 Comparison of Hhlrl Aaiplltude Ratios 

Case 1 Is a typical non-uniform rotor system with Identical plain cylin- 
drical end bearings, L/D ° 1.0, operating with eccentricity ratios 

“ 0.2. The response of this rotor In plain cylindrical bearings 
Is shown In figure 32. To obtain the natural frequency of such rotors 
1n rigid bearings, the following bearing coefficients were used: 

Kyy = Ku = 100 X 10^ Ib/in 

»\A yy 

hy = V “ “-O ’Win 

^yy “ ®»/w ” 0.0 lb sec/in 
AA yy 

“ 0.0 lb sec/in 

The response of the equivalent rotor for this system is shown in figure 33. 

It may be observed that the maximum first-mode stepped shaft response 

-3 

amplitude at mid-span is 2.2 x 10 in and that the mid-span whirl ampli- 

-3 

tude of the equivalent shaft is 5.0 x 10 in. The ratio of these numbers 
is: 


“ = = 2.3, as indicated previously 

\ 2,2 x lO""* 

Case 2 is the same as case 1, but with all dimensions multiplied by 2.0: 
See table 2. The operating eccentricity of both bearings is e = 0.2. The 
rigid bearing response is shown in figure 34, the flexible stepped rotor 
response is shown in figure 35, and the equivalent rotor response in equi- 
valent bearings is shown in figure 36. The ratio of the maximum first 




mode mid-span whirl amplitude of the equivalent rotor to that of the 
stepped shaft is A /A = 20.10”^/8.10''^ = 2.5 in this case. 

Case 3 is the same as case 1, but with all dimensions multiplied by 
3.0. Operating eccentricity ratio is again e = 0.2. Results are shown 
in figures 37, 38 and 39. The amplitude ratio for this case is A„/A^ - 1.7. 

Case 4 shows the influence of two steps on the shaft, with two moderately- 
long cylindr cal shaft sections supported in end bearings. The bear- 
ing eccentricity ratio is again e = 0.2. Results for the rigid bearing 
case, stepped rotor case, and equivalent rotor case, are shown in figures 
40, 41 and 42, respectively. The amplitude ratio is A /A^ = 1.7. 

Case 5 is a double overhung rotor with a flexible shaft between bearings. 

The bearing operating eccentricity is e = 0.2. Results are shown in figures 
43, 44 and 45. The amplitude rat’o for this case is A^/A^ =11.5. This 

result is high because modeling as an end-bearing rotor is unsuited to 

this case. 

Case 6 considers the influence of unsymmetrical shaft sections, and also 
the influence on rotor response of bearings which operate at different eccen- 
tricity ratios. Results are not shown for this and remaining cases in the 
interest of brevity. The amplitude ratio# in this case, was A /A^ =1.7. 

Case 7 has the same rotor as case 6 and the operating eccentricity of the 
left bearing is increased to e = 0.4. The amplitude ratio was A^/A^ = 1.63. 

Case 8 is a modification of the case 6 rotor, with a longer left section in 
the same bearings as case 6. The amplitude ratio is A^/A^ = 1.7. 

Case 9 is the case 8 rotor for which the operating eccentricity of the left 
bearing has been changed to e = 0.4, with the longer left section of the 
case 8 rotor. The amplitude ratio is Ag/A^ = 2.0. 

Case 10 is a three mass rotor with an overhang at each end. The connect- 
ing shaft section is slender. The rotor operates at eccentricity e = 0.2 








for both bearings. The lowest rotor whirl mode is a free-free mode. As 
with case 5, response results for the actual rotor and for the equivalent 
rotor are not similar in form. The amplitude ratio Ag/A = 0.27, pri- 
marily due to the lack of damping in the bearings. This result is not 
regarded as relevant because the end bearing equivalent rotor is too 
strongly affected by bearing stiffness to assume a free-free bending mode. 


5.3 Results 

The main result from these studies was that the equivalent rotor procedure 
appeared to give results for maximum mid-span amplitude ratio A^/A„ which 
lay between 1.63 and 2.50. Thus the equivalent rotor predicted mid-span 
amplitudes which always exceeded those of the stepped rotor by between 
63 percent and 150 percent. The equivalent rotor maximum amplitudes may 
therefore be taken as conservative estimates of maximum whirl amplitude. 

Two studies of overhung rotors were made, cases 5 and 10. In both instances 
the results observed were thought to be unproductive because the overhung 
model did not match the end bearing model. Earlier it was assumed that with 
flexible bearings the free-free mode would dominate the response, and the 
useful comparison would result. It is shown in section 6 that most rotors 
appear to function as overhung rotors, due largely to their overhung 
coupling sections. Care is therefore necessary in interpreting the results 
any rotor in terms of an equivalent end bearing rotor. 

An attempt was made to correlate higher modes and to develop amplitude ratios 
for these modes. The results were not considered sufficiently successful 
to report here and the results are not included. Basically, the modeling 
procedure appears to suffer wherever the system model departs from mid-span 
symmetry. For symmetrical systems only odd-order modes were excited, whereas 
for unsymmetrical systems the even and odd modes were excited. This made 
comparison difficult. 






6. BALANCING STUDY OF SIX ROTORS. 

6 . 1 P rocedure 

Six practical rctor systems for which geometry, operating data and 
unbalance response data is available were examined to provide informa- 
tion on the following questions: 

a) How readily can the response results computed in the present 
study be used to predict practical rotor whirl amplitudes? 

b) How suitable are the flexible rotor balance criteria in 
DP 5343 when applied to specific practical rotor cases? 

c) How suitable are the criteria parameters presently in use for 
rotor balancing, i.e. journal velocity, pedestal velocity? 

Where possible, the practical rotor cases chosen had both test data and 
calculated results data available. Rotor system details are given in 
table 3. Results are compared with existing balance criteria in table 4. 
Additional details used in the calculations reported in these tables were 
obtained from references cited herein unless otherwise stated. 

6.2 Case 1. Lund-Orcutt Rotor. Figure 46. 

Lund and Orcutt [5] have described the calculation and testing of a multi- 
disk rotor in tilting-pad bearings. The correlation shown between their 
test data and their calculated rotor response vs speed data is close. The 
results of the Lund-Orcutt study were subsequently confirmed by Thomas [6] 
The Lund-Orcutt rotor had a single central disk weighing 36.0 lb mounted 
on a uniform shaft weighing 88 lb. The rotor system critical speed is 
given as 12,800 rpm. System dimensions and operating details are quoted 
from ref. [5]. The stiffness parameter C/6 is found as follows: 

1 . Bearing clearance, C 

For preload, m = 0.5 C = 3.75 • 10"^ 

Operating eccentricity at 12,800 rpm: e = 0.2 





TABLE 3. DETAILS OF SIX PRACTICAL ROTORS FROM LITERATURE 



Turbine 















































Table 4. Interpretation of rotor whirl using 150 balance and vibration criteria . 












































2. Mid-span Deflection, iS 

for a rigid-bearing uniform rotor 

For Nj, = 12,800 rpm, 

3. Stiffness Ratio C/6 

C/6 = 9.47 

4. C. G. Eccentricity a = 500 x 10"'^ in (given) 

Figure 22 applies for a uniform rotor in tilting-pad bearings 
operating at journal eccentricity c ~ 0.2. The stiffness curve 
C/6 = 10 may thus be used for comparison with the Lund-Orcutt 

-3 

results for this case. This indicates a whirl radius r = 2.2 x 10 in 

-3 

at the bearings and r = 2,6 x 10 in at mid-span. Values of r/a 
are given below. 

5. Whirl radius ratio (test) 

Bearing r/a = 2.2 x 10"^/0.5 x 10" 

Mid-span r/a = 2.6 x 10"^/0.5 x 10 

6. Whirl radius ratio (calculated) 

Bearing. Figure 22, r/a = 1.6 

Mid-span. Figure 25. r/a = 1.4 


= 4.365 
= 5.158 


6 

6 


5 

3¥4 


Wk 

El 


3.96 • 10‘ 


Calculated values of whirl radius ratio appear somewhat low for the rotor 
configuration involved. This may be due to the stiffness ratio C/6 = 9.74 
which corresponds to a stiff flexible rotor in damped flexible bearings. The 
response charts [5] for this case suggest a more flexible rotor in stiff 
low-damping bearings. Test values of C/6 appear to be correct. Similarity of 
response for the equivalent rotor and for the test rotor may be observed by 




comparing figures 25 and 47. 


A comparison of balancing criteria for the Lund-Orcutt rotor are shown 
in table 4. The Lund-Orcutt single disk rotor is a class 3, category II 
rotor. Rotor category data is given in table 5, from ISO DR 5343, Ref. [11]. 
These criteria values indicate the following conditions for the Lund-Orcutt 
rotor: 

1. Rotor Quality Grade : Small turbine rotor G 6.3 

2. Permissible Vibration Velocity (table 5) 

Journal: V = 0.283 in/sec r.m.s. 

Mid-span: V = 0.567 in/sec r.m.s. 

3. Corresponding Whirl Amplitude at 12,800 rpm 

Journal: r = 0.21 x 10“^ in r.m.s. 

_3 

Mid-span: r = 0.42 x 10 in r.m.s, 

4. Corresponding C. G. Eccentricity (implied by whirl radius) 

Journal: a = r • = 0.21 = 48 x 10 ^ in 

Mid-span: a = 0.42 x = 81 x 10 ® in 

5. Permissible Whirl Radius/Clearance Ratio 

Journal: r/c = 0.21 x 10”^/x 10"^ = 0.11 

Mid-span: r/c = 0.42 x 10"^/x 10 = 0.22 

The required residual C.G. runout implied by the above results is consis- 
tent with practical balance quality for this type of rotor. The importance 
of relating balance criteria to measurement location on the shaft is 
clearly demonstrated by the differing values of whirl radius at the journal 








and at mid-span. It is also apparent that observed whirl values depend 
on the rotor mode shape, Coupling mass causes a rotor to whirl more in 
a free-free mode, whereas mid-span mass tends to cause half-sine modes. 

The implied (whirl radius/clearance) ratio values derived from the DR 5343 
criteria would be consistent with stringent balance criteria, considering 
that they apply to resonant conditions with a sharp peak, figure 47. The 
observed values of r/c in table 5 appear to be more reasonable for this 
condition. 

6.3 Case 2. Rieger-Badgley Rotor. Figure 48, Reference [15], 

This rotor has a non-uniform distribution of mass and stiffness and operates 
in two tilting-pad gas bearings close to its ends. Two critical speeds 
exist within the range tested at 72,200 rpm and at 86,500 rpm. The same 
free- free rotor mode is involved in each instance, and the difference is 
due to bearing properties in the horizontal and vertical directions. 


The stiffness ratio is high because of the high critical speed. This value 

corresponds to a stiff rotor in flexible bearings. This rotor whirls in a 

free-free mode, and this is confirmed by the critical speed ratio (w/w^) of 

2.25 in figure 22 for a corresponding C/6 value of 150. The rotor has an 

unbalance of 500 microinches at each disk, oriented in a planar manner which 

tends to induce the observed free-free whirl mode. Critical whirl radii 
-3 "3 

of 45.0 X 10” in (bearing) and 12.5 x 10 in (mid-span) correspond to 
whirl ratios (r/a) of 90.0 and 25.0 respectively. The calculated (r/a) 
ratios obtained from figures 49 and 50 are 25.0 and 25.0. These values also 
include the effect of higher damping (and stiffness) present in the oil 
film bearings for which the charts were developed. Adjusted values of (r/a) which 
allow for the influence of bearing film compressibility are 47.3 and 6.55, 
respectively. Whirl radius/clearance ratio values corresponding to these 
results are high: (10.5, 1.46) for rotor results and (5.56, 5.56) from 

chart results, at resonance. These high values are due to the high unbalance 
(500 microinches) the prescribed distributions of which couples strongly with 
the free-free mode, as indicated above. 







The balance criteria conditions indicate that t>'*s Is a class 3, category 
IV rotor for which the permissible r.m,s. velocity levels are 0.331 in/sec 
(bearing) and 1.102 in/sec (center). The implied residual c.g, eccentricity 
to achieve these levels based on the corresponding r/a values given in table 
5 are 1.7 microinches and 5.8 microinches. With these c.g. eccentricity 
values the whirl radius/clearance values are 0.03 (bearing) and 0,09 (mid- 
span). 

It is apparent that the DIS 5343 balance criteria appear to be excessive in 
this instance. If resonant r/c values of 0.5 (bearing) was permitted, the c.g. 
eccentricity would be 24 microinches. This is still stringent, but attainable. 


6.4 Case 3. Lund-Tonneson Rotor. Figure 50. Reference [16], 

This rotor system consists of several disks on a flexible shaft, mounted in 
externally-pressurized gas bearings near its ends. The performance of this 
system has been described by Lund and Tonneson [16]. 

The first critical speed occurs at 7850 rpm. The mode shape is not shown 
in ref. [16], but may be expected to be a free- free mode similar to that 
of the Rieger-Badgley rotor. The corresponding system ratio is 2.26, 
which indicates that the rotor is flexible and operates in rigid bearings. 

The stated residual unbalance is 0.095 oz. in, or a c.g. eccentricity of 
161.7 microinches for the rotor weight of 36.7 lb. At the critical speed 
the whirl radius is 0.36 x 10 in at the bearing, and 1.0 x 10“ in at 
mid-span. These figures correspond to whirl ratios (r/a) of 2.24 at the 
bearing and 6.18 at mid-span. For comparison, plain cylindrical bearings 
operating at eccentricity e = 0.2 with an equivalent rotor C/6 =3.0 the 
whirl ratio is 1.30 at the bearing and 13.0 at mid-span. 

The test rotor resonant r/c values are 0.23 (bearing) and 0.64 (mid-span). 
The equivalent rotor values are 0.13 (bearing) and 1.34 (mid-span). Both 
conditions appear to correspond to a medium-unbalanced condition which 
the 161.7 microinch eccentricity would suggest. 




The balance criteria in table 5 indicate that this is an ISO class 3 
category II rotor for which the permissible r.m.s. velocities would be 
0.47 (bearing) and 0.89 (mid-span). The corresponding r.m.s. whirl 
amplitudes at 7850 rpm would be 0.58 x 10“^ in (bearing) and 1.08 x 10"^ in 
(mid-span). Using the whirl ratio, these amplitudes imply C.G. eccentri- 
cities of 366.2 microinches (bearing) and 247.1 microinches (mid-span), 
respectively. The corresponding whirl radi us/clearance values are 0.52 
(bearing) and 0.97 (mid-span). 

The results show that the ISO DR 5343 criteria values are fully met with 
i;he rotor balance stated in table 5. For gas bearing rotors stringent 
balance control is necessary, which may account for the low residual 
unbalance observed. 


6.5 Case 4. 6as Turbine Rotor. Figure 52. 

This rotor weighs 32.0 pounds and is mounted in fluid-film tilting pad 
bearings. The first bending critical speed occurs at 7500 rpm though 
the operating speed range extends through 52000 rpm, and includes several 
higher critical speeds. The mode shape is again a free-free mode. 

The ,'..g, eccentricity of 500 micro inches causes a critical whirl radius 
0.2 X 10 in at the bearing and 0.5 x 10 in. at midspan. These values 
correspond to whirl ratios r/a of 0.4 (bearing) and 1.0 (bearing). The 
corresponding clearance ratios r/c are 0.125 (bearings) and 0.313 (midspan). 

The stiffness ratio for this rotor system^based on the bending critical 
speed of 7500 rpm> is 2.03. This indicates a flexible rotor in stiff 
bearings. From figures 22 and 25 for til ting-pad bearings the response 
of the equivalent rotor at resonance is r/a = 1.4 (bearings) and 3.9 
(midspan). The corresponding clearance ratios r/c are 0,44 (bearings) 
and 1.22 (midspan). 

The balance criteria in table 5 indicate that this is a Class 3, Category 
IV rotor for which the permissible rms velocities are 0.33 in/sec (bearing) 




and 1.10 in/sec (midspan). The corresponding rnis critical speed amplitudes 
are 0.42 x lO’^ in (hearing) and 1,4 x 10"^ in (midspan). The implied c.g, 
eccentricity which would cause these whirl amplitudes are 1485 micro inches 
(bearing) and 2008 (midspan). 

From this is appears that the existing balance condition of 500 micro- 
inches is adequate, both from existing criteria and from the clearance 
whirl ratio values, 


6 . 6 Case 5. Kendig Gas Turbine Rotor. Figure 53. [..fer e nce [17 ], 

This is a small gas turbine rotor weighing 50 pounds and operating in 
til ting-pad gas bearings. The service speed is 66,000 rpm. There is a 
damped rigid body composite critical speed at 10,000 rpm and a bending 
critical speed at 24,950 rpm. Modes for both these critical speeds are 
shown in figure 54. The stiffness ratio c/6 = 10.5 is based on the 
bending critical speed and indicates a flexible rotor in flexible bearings 
at this speed equivalent rotor response chart for this stiffness ratio 
shows a damped critical speed peak at approximately 40 percent of the 
free-free critical speed. This result supports the equivalent rotor 
model for this case. 

Calculations were made to determine the response of this rotor system to 
a residual inbalance of 0.3 oz. in or 375 microinches c.g. eccentricity. 
Whirl radius of 7.0 x 10" in (compressor bearing) and 9.0 x 10 in 
(turbine bearing) were indicated at the bending critical speed, and 
lower whirl amplitudes at the rigid body critical speed. These values 
correspond to whirl radius ratio values of 24.0 and 12.0 respectively. 

The equivalent rotor whirl ratio values are 5.3 (compressor) and 5.3 
(turbine). The corresponding critical speed clearance ratios r/c are 
12.0 and 6.0 for the calculated rotor and 26.5 and 26.5 for the 
equivalent rotor. 

The criteria table shows this rotor to be a Class 3 Category IV rotor for 
which the permissible rms velocities are 0.331 in/sec (journal) and 




1,10^ in/sec (intcispan). Pamissible cHffeal whirl aiiiplltudss are 0.26 x 10" 
in/sec and 0.90 x 10"^ in/see. These amplitudes imply c.g. eccentricities 
of 15,3 microinches and 106 niicroinches. The corresponding clearance 
ratios are 0,49 and 1.7. 

Practical considerations require clearance ratios of 0.5 at the bearings. 

This corresponds to a c.g. eccentricity of 15 niicroinches. This fine 
balance is consistent with operation through a bending critical speed with 
gas- bearing turbomachinery. 


6.7 Case 6. Kendig Steam Turbine Rotor. Figure 55. Reference [17 1. 

This rotor system consists of a 5000 1b. medium power steam turbine 
rotor in end bearings. Critical speeds exist at 2300 rpm (rigid body) 
and at 12,000 rpm (free-free bending). The stiffness ratio c/6 = 14.5 is 
based on the bending critical speed. Reference to figure 22 for an 
equivalent rotor in tilting pad bearings shows a rigid body critical 
speed around w/w^ - 0.5 and a free- free mode at 2.25. 




7. COMMENTS ON FLEXIBLE ROTOR BALA.CINC CRITERIA 


Details of the specific rotors studies are given in Table 3. Details of the com- 
parison with ISO Draft flexible rotor balance criteria are given in Table 4. The 
criteria numbers are obtained from the ISO documents indicated, All rotors are 
flexible (Class 3) rotors. The ISO rotor category refers to the type of machine 
in which each rotor is used, as indicated in Table 4 of ISO DR 5343. 

The balance criteria values given in ISO DR 5343 are based on VDI criteria, which 
are compatible with criteria for vibration in rotating machinery given in ISO 
DR 2372. The permissible journal rms velocity values in Table 4 apply to each 
rotor. Values are given for the journal, and at midspan using the correction 
factors in Table 4 of ISO 5343. Permissible amplitude values were obtained by 
using v - wr where v is the rms velocity, w the rotor speed in rad/sec and r 
is the permissible rms amplitude. Using the uniform rotor response charts 
allowed the implied c.g. eccentricity to be deduced corresponding permissible 
values of (whirl radi us/bearing clearance) were then obtained for each bearing. 

Table 4 shows that using the ISO criteria with the response charts gives 
permissible resonant whirl radius ratios (r/c) between 0.03 and 0.72 at the 
journals of these six rotors. The exception to this is case 5. The small gas 
turbine rotor, which shows excessive values at both its critical speeds. This 
discrepancy is due to excessive rotor unbalance. Good balancing could reduce 
the c.g. eccentricity to 25 x 10" in, which would in turn reduce the r/c values 
to 0.6 and 2.0 respectively. The latter value is still large, possibly due to 
inaccuracy in defining the peak accurately in the calculations, 

Table 3 shows results for critical whirl radius r at the journals and at midspan 
obtained from the results of several authors, under critical speed conditions. 
These values range from 0.2 in to 70.0, depending on the amount of residual 
unbalance present. Normalizing these results with c.g. eccentricity, a, gives 
(r/a) values between 0.4 and 24.0 (journal) and 1.0 and 167 (midspan). Comparable 
results obtained from the equivalent rotor calculations show values between 1.3 
to 25.0 (journal) and 1.4 to 2b. 0 (midspan). In general, the rotor r/a values 
may be up to four times larger at the bearings than the values of r/a determined 
by the equivalent rotor method. As these are resonant r/a results, it is 
possible (perhaps likely) that peak r/a values in practice would be smaller than 
indicated in the actual rotor results, due to non-linear effects in the bearings. 
At midspan, the correlation is again similar with the actual r/a values 






exceeding the equivalent rotor r/a values by ab^i 4:1. 

During the analysis of results It was noted that most actual rotors behaved more 
like free- free beams In their bending modes than like pinned-pinned beams . This 
Is apparent when It Is observed that most rotors have overhung end drive sections 
carryl ng couplings. It Is felt that If the equivalent rotor method were based 
on results from overhung rotors, closer correlation could have been obtained In 
all results. The extension of the equivalent rotor method to Include a free-free 
beam analysts is recommended for further study. 

The accuracy of results obtained in this study was also affected by the impreciseness 
of the resonant amplitude values both In the actual rotors and ai the equivalent 
rotor calculations, As the final results are strongly dependant on both peak values 
it is evident that discrepancies and some randomness will occur in the correlation 
of results. As a present guideline, If the equivalent rotor method described 
herein is used, and the results for r/a are multiplied by a factor of 4, this 
appears to allow sufficiently consistent resonant r/a values to be defined for use 
in practice. When these values are then utilized with the r/c values at the bearings, 
a decision on the quality of balance required can be made. 

The final conclusion is that the proposal equivalent rotor method appears to 
offer good promise as a simple procedure for deciding the quality of balance 
required in a given rotor, but that the procedure should more likely be based on 
the lowest free-free mode of a uniform flexible rotor in bearings rather than on 
end-bearing rotor modes, as was done herein. When this is done the method 
could represent a consistent basis for balance quality selection. 






8. CONCLUSIONS 

1. The equivalent rotor method represents a practical procedure for selection 
of flexible rotor balance criteria, based on the uniform rotor/specified 
bearing charts presented in this report. 

2. General unbalance response charts for uniform rotors in specified fluid- 
film end bearings have been developed. Data is available for rotor 
response at mid-span and at journals for each bearing type. These 
charts are suitable for specifying balance criteria for flexible end- 
bearing rotors. 

3. Use of these general response charts should extend to end-bearing rotors 
only. Application of the charts to rotors with end overhangs may give 
results which are in error. 

4. The charts apply best to end-bearing rotors with mid-span symmetry. 

Where system symmetry decreases for any reason, the correlation between 
actual response amplitude and equivalent rotor response amplitude also 
decreases. 

5. Difficulty arises in specifying general balance criteria for flexible 
rotors which operate through more than one critical speed because one 
mode may not be excited as strongly as another mode by the rotor 
residual unbalance distrlfcjtion. 





9. RECOMMENDATIONS 

1. Similar response charts should be developed for a uniform rotor with 
overhangs. This would provide general response data for the two 
main classes of flexible rotor systems. 

2. Practical data should be obtained from flexible rotors operating in 
the field unbalance to demonstrate the suitability of specifying 
balance criteria based on the equivalent rotor procedure. 
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n. APPENDIX - COMPUTER PROGRAM ROTOR 


11.1 Program Capabilities and Limitations 

ROTOR Ig a genoral purpoce computer program for 

the unbalance reoponse analycis of a uniform elastic 

rotor supported in fluid ~film bearings. The program 

is capable of assembling a rotor with any or all of 

the following basic components, 

1, Beam elements- having distributed mass and elastic 
properties and constant E,A.and 1 throughout the 
length, 

2, Disks- with concentrated mass MO, polar and 
transverse mass moments of inertia Ip and 
respectively, 

3. Fluid-film bearings- described by eight speed 
dependent dynamic stiffness and damping properties 

^XX* ^Xy» • • • ♦ D D . . . , 

4. Unbalance forces- represented by the real and 
imaginary components of a rotating force vector 
at each node. 

The maximum number of elements are 7 thus fixing 
the maximum number of nodes at 8, At each speed 
increment the complex nodal displacements and 
rotations are calculated and used in formulating the 
major and minor ellipse radii and the ellipse angle 
from the positive x axis to the major axis of the ellipse 
and the program gives this information as output. 










Limitationo of the program are i 

1, Only two identical bearings may be accurately 
represented. For multibearing use, the 
statically indeterminate support problem must be 
solved first. 

2, Short stubby sections heavily loaded in shear 
may not be accuraf'ely represented. 

3, The same weight density and Young’s Modulus is 
used for each beam section, 

4-, Static deflection due to gravity is not taken 
into concideration thus limiting the analysis 
to the classical "vertical rotor" problem. 

5. Steady state response is only concidered. 






n.2 General Programing Information 

ROI'OR ie written in Fortran IV and wao developed using 
a Xerox Sigma 6 computer. Data input is from a card reader 
whoso device unit no, is (105) ^ind output is to a line 
printer whose device unit no. is (108), All computations 
are executed using complex double precision arithmatic. 

The main program plus eleven subroutines contained ^?5 
Fortran statements. The program requires 15.8 K words 
of main computer core. 

The functions of the main program and of the subroutines 
are as follows, 

Ip .main program - array declarations, reads in and writes 
out data input for checking, acts as an executive routine 
for calling subroutines, assembles the structural stiffness 
matrix, generates the force vector and solves for the 
major and minor whirl ellipse axis. 

2, ELEM - sets up an 8x8 dynamic stiffness matrix for a 
single beam element, all terms in the matrix are real 
double precision. 

3, EDISlCCi - sets up an 8x8 complex double precision matrix 
which reflects the effects of the addition of a disk 

to the left end of the element. 

EDISKR - same as 3 except for a disk added .j the right 
end of the element. 
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k 

f 




5, EBBAIili - seta up an 8x8 complex double precioion matrix 
which reflects the effects of a fluid-film bearing added 
to the left end of the element, 

6, EBEARR - same as 5 except for a fluid film bearing 
added to the right end of the element. 

7, ZEROM - initializes a real matrix to zero. 

8o ZEROMC « initializes a complex matrix to zero, 

9, CADDMl - adds a real matrix to a complex matrix. 

10. CADDM2 - adds two complex matrices. 

11. GMULTM - multiplies .two complex matrices 

12. CINV - inverts a complex double precision matrix 
using the Gauss Elimination with partial pivoting. 










n.3 Input Data Format 


Data input to ROTOR is in tho form of punched cards. Seven 
sets of data (1-7) are required, with the number of input 
cards per set depending on the particular problem being 
solved. The definition of the input parameters, the order in 
which they should appear and their format is as follows. 

DATA SET 1 - General element information 
One card (2I5,3F20,3) 

NEliEM - total no, of elements used (max. 7) 

NODES - no, of nodes ' (max. 8) 

EMOD - Young's Modulus Ib/in^ 

ERHO - weight desity Ib/in^ 

DATA SET 2 - Speed information 
One card (3P20,3) 

BSPEED - begining speed rpm. 

SPEEDi’- speed increment rpm. 

PSPEED - final speed rpm. 

DATA SET 3 - Unbalance forces ' 

One card for each node (3F20.3) 

CUBALP - cos component of unbalance oz-in. 

SUBALF - sin component of unbalance oz-in, 

DATA SET ^ - Control cards and specific element information 
Two cards for each element, a total of 2xN£LEM cards 
card 1 ( 515 ) 

lELEM - 1 if element exists 

0 if element does not exist 







IDISKL - 1 if loft end diok ie present 
0 if no left end disk 
IDISKR - 1 if right end disk is present 
0 if no right end disk 
IBEARL - 1 if left end bearing is present 
0 if no left end bearing 
IBEARR - 1 if right end bearing is present 
0 if no right end bearing 
card 2 (2I5,3F20.3) 

NA - left end node ‘.of element being specified 
NB - right end node of element being specified 
EINER - element inertia in^, 

EAREA - element area in^ . 

ELEN - element lenght in, 

DATA SET 5 - Disk information 

One card for each node, (3F20.3) 

DW - disk weight lb. 

DRAD - disk radius in. 

DLEN “ disk lenght in. 

DATA SET 6 - Bearing locations 
One card (515) 

NBL - left end bearing node 
NBR - right end bearing node 
DATA SET 7 - Bearing properties 

Two cards for each speed increment (4F20.3) 








T * 


card 1 

SXX « bearing Gtiffnese in load direotion Ib/in, 

SXY “ crooG coupled stifCneoc Ib/in* 

SYX - cross coupled stiffness 

SYY « bearing stiffness perpendicular to load direction Ib/in. 
card 2 

DXX ” bearing damping in load direction lb. sec, /in. 

DXY - cross coupled damping lb. sec. /in, 

DYX ~ cross coupled damping lb. see, /in. 

DYY - bearing damping perpendicular to load direction lb see/in 






11.4 Sample Input-Output Problem for P.0T0R 
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Af'ALYSrrj 


c 

c 

c 


Cn.visii'Tt; R PPn(jRA.M PHTOR , * . , . 
p-n» TI-? SfKADY STAfF UNPAIiAMCF: RtJ.f^Prt'.rF 
,'jF p-Tr-Bg ri Ff.utD-PiLF csFAPrias. i*p w nf Royop 
f*:rP.wF*'T« *«AV mF: axially AvSSF^.H{ &:d. THP QHT'W 

f.;i F^t.?.',' *„Agi.S a - n ELAS'fir PpnpFkTIp;? APF riSTwlRrifFn 
:;\|B‘tyWf rY ar nfifi THP; LK^jGTH Of 'jl-’t. KLpA'HlI . 

L’vpai.AniCr, FOPCKS, FLUlD“rrLP »aa«IfLS AM|: DISK.H 
ARP Air, nWF.R AT EACH NODAL POlMf. AT F.ACH SPEED THE 
OYHayiC STIFFHESS MATRIX IS FORMED AND INVERTED BY 

GAUSS elimiwation with partial pivoting, at each 

NOPAL POINT THE STEADY STATE WHIRL ORilt IS 
CALOULATEO DUE TO SPECIFIED UNBALANCE IN THE SYSTEM, 

THE unbalance FORCES ARE BEPRFSEHTED RY A ROTATING 
VECTOR. THF FLUIO*FILM BEARINGS ARE REPRESENTED BY 
EIGHT SPEED DEPENDENT DYNAMIC STIFFNESS AND DAMPING 
COEFFICIENTS. THE DISKS POSSESS MASS AND GYROSCOPIC 
COUPLING OPOPERTIES. COMPUTATIONS ARE PERFORMED IN 
.. complex double PRECISION ARITHMATIC, 

0 ^ ^ ^ ;<o!s ^ ^ ^ ^ ^ ^ ^ sSi ){«!? 9 ¥ J|! ^ ^ ^ K! ^ iK 5(« 

C C.B. THOMAS, MjF.^RIEGER, A. ARvSLANCAN 
C wnCHFSTFR INSTITUTE OF TECHNOLOGY JUNE, 1974 

C MAIN PROGRAM 

C ^ :j£ ^ ^ :t« ^ # 5*< Sj! ^ 4! >i= ^ ^ ^ 9 S>t 2^ ^ !!< * »K >J{ ^ J}! 4i 

I*-'PLICIT REALMS CA’’H,0-Z) 

PFfiL=f-‘8 EELKC9 .8) ,CUBALFC8) ,SUBALFC82,CENTFCCS) ,CENTFS(83 
REALM8 DRADC85 ,DLENC8) ,DWC^) ,EINERf7 5 ,EAREAa3 Er.FWC75 
rOMPLEX^lS EDKLC8,8) ,EDKRC8 .8) .FBKLC8!fi) ,EBKRC8,8) 

COMPLEX’!* Ifi 5DKC3P,82) ,SDKlC;h,32? ,FVE6( 3^ 1 3 ,XX oL 1) 
DIM|N|ION IELEMC75gIDISKLC7),IOISKRC73 ,IBEARLC7),IBEARR(7) 

CGMFriN EMOD . E XNEP , EAREA , ELEN , ERHO 
common f elk , EDKL , tDKR , E8KL , EBKR 
common DW,DRAD.DLEN 

common SXXL,SXYL,SYXL,SYYLpDXXL,DXYL,DYXL,DYYL 
COMmhn SXXR,SXYR,3YXR,SYYR,DXXR,DXYR,DYXP,DYYR 
c if: i«« ^ ^ ii< :<< i|' ift ^ >)t ){t )(s Sjc :j< ^ J|( s;! St M ^ sx i(< M 5{< M « si* ^ ^ 

C READ IN AND WRITE OUT INPUT DATA 

(j 1|{ !j! Sft !!s )|{ ){< lit ^ 5{« # }}! >|t ^ ^ )jt :{t « J|t ^ ^ J{! !,t *)}! ^t ;){ M ^ « s’.* ^ ^ :;< * !^t * 

READClOb,!) NELEM, NODES , EMOD , ERHO 
READ C 1 OS. 2 3 BSPEEO , SPEED! , FSPEED 
“ ■ .NODES 

CUBALFCI) ,SUBALF(I3 
, NELEM 

1ELEMCI3 ,IDISKL(I3 .IDISKRCI3 .IBEARLCI) .IBEARRCD 
NACI3 .NBCI3 .EINEPCI) .EAREACI3 ,ELEN(I3 
, NODES 

DWCI3 ,DRAD(I),DLENa3 
N 0 L , N B R 


1000 


1050 

2000 


10 


DO 1000 1=1 
RF,ADri05,23 
nn 1050 I=l 
READ Cl 05, 3 3 
READCIOS, 13 
OQ 2000 I=l 
READa05,2) 

READ a 05, 3 3 
NSIZE=4#CMELEM+13 
WRITEC108, 103 


FORMATC//.1X, 'M0„ ELEMENTS' ,5X, 'NO. NODES SX YOUNGS mu 
15X, 'WEIGHT DENSITY' ,5X, 'FIRST SPEED ', 5 X ,' SPEe6 INCREMENT 
2SX, 'FINAL SPEED' 3 

WRITE(108,153 NELEM, NODES, EMOD, ERHO, BSPEED.SPEEDI, FSPEED 
15 F0RMATC/,2X,I5, 10X.I5,8X,D11.5,9X,F5,3,12X,F7.1, lOX, 
12X#F7«1) 

WRIT S' (108, 20 3 

20 FOP.V.ATC//, IX, 'ELEMENT NO . ' . 5X , ' NODE A',5X,'NODE B',5X, 

1 'ELFMENT INERTIA' ,5X, 'ELEMENT AREA' ,5X, 'ELEMENT LEN(5tA ' 3 
DO 2050 1=1, NELEM 

2 0 50 WRITE Cl 08, 25 3 I . N A C 1 3 , MB (1 3 , FINER ( 1 3 , EAREA C 1 3 , ELEN (1 3 
25 FOPMATC/.5X,I2,i2X,I2,9X,I2 9X,D10.4,9X,D10,4,9X,D10.43 
WRITEC 108 , 30) 

30 PORMATC//, IX, ' UHALANCE FORCE ', 5X NODE NO.' 
irn.e; tompomitmt i . c;y . i siN COMPONENT '3 


MOD' 


1 'COS COMPONENT' ,5X, 

DC 3000 1=1, NODES 
3000 '‘/PITFC 108 . 35 3 I , CU8ALF ( 1 3 , SU6ALF ( 1 3 
35 F0PMATC/,23X,I2, 12X,F7,3, 12X,F7.33 
WRITFC108, 37) 

37 FORMATC//, IX, 'DISK INFORMATION ', 3X 


5X 


'NODE NO. ' ,5X 


/ ^VA#'UJLV. 

I'DISK WEIGHT' ,5X, 'DISK RAD. ', 5X ,' DISK LENGTH') 
DO 3010 1=1 .NODES 

ORIGINAL PAGL IS 

OF Pl^OR QUAJ ,!*F. 




3010 '.vRTTG-(l08.3»?) TiOaCD.DRADCI) ,0I.PN(X) 
38 rnp .iATC/,i3X,I2ilOX,F8.3,6X,Ffi. ■ “ - ~ 


HOOej ' ,I5,SX 


4,7X,F8.4) 

WRTIRa08,40) NRLf-iaR 
40 P’ORwarC//, IX, ' LEPT SND RFAPIMU 
I'PinHT EMO OKARI^JG ^0DF:',TS) 

r PFGT'V^'iMf, cF SPP.eD r.onp 

444 POVFGAs8SPPEO*2 .OOOI'B. 14lS<}26bnn/60,0C/0 
READ n 05, 4 ^ SXXL,SXYL,SYXL,SYYL 
READn05,4^ DXXL,DXYL,DYXL,DYYf, 
RFADC105,4) SXXR,SXYP,SYXR,SYYR 
PEADCi05.4) DXXR.OXYR,DYXR,DYYR 
V!RTTF( 108,45) 8SPSED 
FOPMAT(///ilX, 'RniDR SPEEDS' ,FP 


45 


1) 


, 'KPL=' ,F10, 
,F10.n 


JC///.1X, . 

WRTTFC 108,50) SXXL,SXYL,SYXD,5YVL 

50 FQR^.AT(//,^X, 'DEFT BEARING INFDRMATION ' , 5X 
1 'KXYL=' ,Fina,5X, 'KYXLs' ,F10.l,5X, 'KYYLs' 

WRtTEC108,55) nXXL,DXYL,0YXL,DYY j 
55 PORMATC/,30X, 'DXXLs' ,n6.1,5X, ' nXYLs ' , FIO . 1 , 5X , ' DYXLs ' , 
1 SX, 'DYYLs' .FlO.n , 

WRI r EC 108,51) SXXR,SXYR,SYXR,SYYR 

51 FORvarC/, IX. 'RIGHT bearing INFORM aTION 4p ' KXXR= Fl 0 , 
1 'XXYRs' ,F10.1,5X, 'KYXRs' ,P10.l,5X, 'KYYRs' ,F10.1) 

WRITEC108.52) DXXR, OXYR ,DYXR , DYYR 

52 FORMATC/, 30X, 'DXXR=' , FIO , 1 , SX , ' DXYR= ' , F 10 , 1 , 5X , 'DYXR=' , 
1 5X, ' DYYRs' ,F10,l) 

C INDIVTDUAD ELEMENT SET UP FOR DYNAMIC STIFFNESS 
C MATRIX USING INPUT CONTROL 

DO 3050 Jsl.NELEM 
IFCIFLEMCJ)) 200,200,100 
100 CONTIiNtlE 


CALL ELEMCJ,R0MEGA) 

GO TO 2Q1 

200 CONTINUE 

CALL ZEROMCEELK.8,8) 

201 irCiniSKLCJ) ) 400,400,300 
300 continue 

CALL EDISKLCJ.ROMEGA) 

GO TO 401 

400 CONTINUE 

CAIL ZEPOWC(EDKL,8,8) 

401 IFCIDISKRCJ)) 600,600,500 
500 CONTINUE 

CALL EDISKRCJ,ROMEGA) 

GO f(2 601 

600 continue 

CALL ZEPOMCCEDKR,8,8) 

601 TFCIREAPLCJ)) 800,800,700 

700 continue 

CALL EBEARLCROMEGA) 

GO TO 801 

800 CONTINUE 

CALL ZER0MCCEBKL,8,8) 

801 IF(IBEAPPCJ) ) 900,900,850 
850 CONTINUE 

call EBFAPPCROMEGA) 

GO TO 901 

900 CONTINUE 

CALL ZER0MCCE3KR,8 ,8) 

901 CONTINMIE 

CALL CACDUl CEELK,E0KL,8,3) 

CALL CAr;DM2CE0KL,EDK .8,8) 

CALL CACDM2CEDKR,£Hf<i., ,8,8) 

CALL CACD*"2(FBKL,EBKP,8,8) 

TCCJ.GT.n GO TO 902 

C initialize TME structural dynamic stiffness MATRIX 

C TH 7,KR0 A‘'D FIT IN THE INDIVIDUAL ELEMENT MATRICES 
CAT.L ZERO’'CCSDK, 32,32 ) 




l,5X, 
FlO.l, 
1 ,5X, 
FlO.l , 


0 02 


<1050 




4000 

3050 


Tisn 
0J = T1 

on 4000 1^ = 1, R 
?-ir. 4050 Dsl.R 

«t‘KaiiOJ3=i.SDK(II,vJJ)+FBKRCK 

•JU = T1 
TT = .[T4-1 
COor r\<IF 

CALL ZF.PnMCCFVEC,32, n 

C BUIIiD THE PARTICULAR i^QPCE VECTOR ASSOCIATED 
C WITH THF problem — ” 

DO 5000 0=1, MODES 

CENTFCCJ)=CUBALECJ)*POMEGA*>ts2/(306.4DO=)‘l6.ODO) 
CEMTFS C J)=SURALFC J)=««R0MEGA'^*2/ C 386. 400^16. ODO) 
B'VECC4*J"3,13=:DCMPLXCCENTFSCJ) ,-CENTFCCJ) ) 
soon B’VKC(4*i'J"l,l)=DCHPLXCCENTFC(J3 ,CENTFS(J3 ) 


e 

C 

C 


IMVERSIOty CF THE STRUCTURAL 
matrix Bi'iD , SOLUTION FOR THE 
AMD R'^TATICMS 

call CTMVCMSIZE/SDK ,SDKI) 


DYNAMIC STIFFNESS 
NODAL DISPLACEMENTS 


CALL CMULTMCSDKI ,NSIZE,NSIZE#FVEC, 1 ,XX3 


TO 880 


C CALCULATION OF THF MAJOR AND MINOR AXIS OF THE 
C WHIRL ORBIT and ANGLE ALPHA FROM THE POSITIVE 
C X AXIS TO THE MAJOR AXIS OF THE ELLIPSE IN THE 
C OIRECTIOm of ROTATION 

no 5050 K=l, NODES 
A = DREALCXXC4*K-3,1) ) 

R=DIMAGCXXC4*K-3,1) ) 

F=DPFAL(XXC4>i‘K-l , 1) 3 
F=QImaG(XXC 4*K-1 , 1) ) 

A1 = A^‘*2 + B=!<*2 + E**2 + F^^2 
A2 = Al*’*<2 

A3 = ~4 .0>i'f-A'«‘F + B*E3 **2 
IFC CA2+A.13 .LT.0.0) GO 
A4 = 0.5>«nSQRTCA2-!*A3) 

GO TO 881 
A4=0,0 

AM = cSQRTC0.5=f‘Al+A43 
AMM = DSORTC0.5!«<A1-A4) 

AI..PHAs: 57,2957E0«ATANC2.0E0#CA*E+B*F) /CF,*E+F=('F 
WRITE Cl 08 ,60) 

P-ORMATC// , 1 X, ' WHIRL ORR IT ' , 5X , ' NODE NQ.', 5 X, 
ISX, 'MINOR ELLIPSE R AD 5 X ELLIPSE ANGLE' 
'WRITE ( 108 ,65) K , AM , AMM , ALPHA 
FQRMAT(/,19X,I2,10X,D15,8,9X,D15.8 ,9X,F8.3) 

continue 

'3SPEED = PSPEED+SPEEDI 
IFCHSPEED.GT.FSPEED) GO TO 9999 
GO TO 444 

C END OF .'^PEED LOOP 


B80 
88 1 




60 


' MAJOR 


65 

5050 


1 
2 

3 

4 

9999 


FORMAT(2T5,3F20 

FnRMATC3F20.3) 

'J’ORMATCSISI 

FQRMATC4F20,3) 

STOP 

END 


3) 






R:«B) )/2.0E0 
ELLIPSE RAD 


EuRM.... SETS UP THE HEAL DnUPLE PRECISION 
C STIFFNESS MATRIX FOR AN ELEMENT 


C SUHRDI'TINE I 

C yXH DYNA.v'IC OX L r r •’* C#0 0 » H l L a r v • rs ►M* C. U QI” C. ^ 

SURROUTTNE ELKMCvI,PnMEGA) 

T.-'PLrCIT HEALER CA-H,n-Z) 

C PAR DYNAMIC .STIFFNESS MATRIX 8X8 REAL 

REALMS EIi''ERC7I ,EAREA(7) ,ELKNC7) 

RFAL>f‘B EFLKC8,P) .DWC81 iDRADCPJ ,ULEri(8) _ „ 

COMPLEX’)' 16 EDKL(8,P) ,EnKR(9,8) ,FRKL(8,8) ,ERKPCR,8) 
common FM(JD,EINER,F.AREAj,ELEM,fc,RHO 
common £ELK,EDKL.EDKR rERKL,ERKP 
common DW.DRAD.DLEN 

common SXXL,SXYL,SYXL,SYYL,DXXL,DXYL,DYXL,DYYL 
COMMON SXXR,SXYR.SYXR,SYYR,DXXR,DXYR,DYXR,DYYR 
ELAM4=CERH0»EAREA(J)’)‘(R0MEGA*’I'25 ^/CEMOD^EINERCO’XSB^ 
ELAM2=DSQRT(ELAM4) 

ELAH=DS0RTCELAM21 

ELL=ELAMX«ELEN(J) 

SELL=DS1NCELL) 

CELL=DCQSCELL1 

SHELLsDSiNHCELL) 

CHFLT, = DCOSH(ELL) 

F1 = SFL.L’(‘SHELL 

F3=CF.LL*CHELL-1 .ODO 

F5=CELL*SHELL"SELL*CHELL 

F6=CELL*SHELL+SELL’i‘CHELL 

F7=SELL+SHELL 

c'8 = SELL-SHELL 

F10=CELL-CHELL 

EC0M = EMQDMEINER(J)*ELAM!)‘*2/F3 
CALL ZEROM(EELK,b\ 8 ) 

EELKC 1 , n=ECOM’)«C-ELAM’l<F6) 

FELKa,2)“ECOM:*(-Fn 
EELKCl ,5)=ECOM’«<(ELAM>l'F7) 

PELKC 1 ,6)=ECOM*(F10) 

FELKf2, n=EELKCl ,2) 

FELK C2 ,2) = ECOM>S<CF5/ELAM) 

FELKC2,5)=-EELKC1,6) 

F.ELKC2, 61=ECOM’)'(F8/ELAM) 

EF,LKC3,3)=EELKCl,n 

EELKC3,41=EELKC1,2) 

EELKC3 ,7)=EELKC 1 ,5) 

EELKC3,8D=EELK(1,6) 

FELK C4, 3)=EELKC2, t ) 

EELK(4,4)=EELKC2,2) 

EELK(4,7)=EELK(2.5) 

feLK(4.8)=EELKC2,6) 

FELKC5, n=EELK(l ,51 
EELKC5,21=:EELKC2,5) 

EELK(5,5)=EELKC3,31 

EELKC5.61=-EELKC3,4) 

FELKC6, n=EELKCl,61 
EELKC6,2)=EELXC2,61 
EELKC6,51=EELK(5,6) 

E£LKC6,6)=PJELKC4,41 

EELKC7,31=EELKn,7) 

FELK(7,4)=EELKC4,71 
fflKC7,7):.”EELKC5,51 
EELKC7 ,81=EELKC5„S1 
FfTLi<(e,31=EELKC3,8) 

FELKC8 ,41=EELKC4,81 
FELK(8,7)=EELKC7,81 
EELKC8,8)=:EELKC6,6) 

RETURN 

END 






,OD01 






c l.FFT Ffh 01 ^K FFK'ECTS 8Xrf COf-PLEX DOUBLE Dt^ECKBTQN 

.SUBfciltiul EDISKLCJ ,ROMfGA) 

T'-PLTCIT PRAL #8 CA-H.O-Z) 

, PEAL^g EJNERC 7 ) ,EAREAC 7 ) ,CLRMC 7 ) 

PEAL*8 EELKC8,8).DWC8) ,DRAn(8) ,DL&Ni - 
COMpLEX'i'U EDKLC8,R) .EDKRC8,8) , EBKL C 8 , 8 ) , EBKR C 8 , 8 ! 

COMMON emod,eimer,earea,eden,erho 

COMMON EELf< . EDKL , EDKR , EBKL , EBKR 
common DW,DRAD,DLEN 

common sxxl,sxyl,syxl,svyl,dxxl,dxyl,dyxl,dyyl 

common SXXR,SXYR,SYXR,SyYR,DXXR,DXYR,DYXR,DYYR 
DMASv 8 = DWfa)/ 386 . 4 D 0 

TINEP“CDMASS/12.0D0)*(3.0D0*DRADCJ)*=<'2 + DLENCiJ)**2' 
PINER=CDMASS/2.0D0)«CDRADCJ)’»‘*2) 

CALL ZERQMCCEDKL, 8 .8) 

fDKLC 1 , l)=DCMPLX(-DMASSif<ROMEGA>«<=t'2,0.0) 

>^DFLC3,3)=EDKLa,n 

EDKLr 2 ,? 3 = 0 C-MPLXC-TIMER*R 0 MEGA=to)c 2 , 0 . 0 ) 
EDKLC 4 , 4 )=EDKLC 2 , 2 ) 

FDML(2,4)=DCMPLXCO.O,-PINER»!ROMEGA>!'=»<2) 

Pr)KLC4,2)=DCMPLXCO.O,PXMER^ROMEGA>f<^2) 

RETURN 

END 


C RTGMT END DISK EFFECTS 8X8 COMPLEX DOUBLE PRECP3SI0N 

SUBROUTINE EDISKRCJ.ROMEGA) 
implicit REALMS (A-H.O-Z) 

REALMS EIMEPC7) .EAREAC7) ,ELEN(7) 

PEAL* 8 EELKC8,8) ,DWC8) ,DRAD(8) ,DLEN(8) 

COMPLEX*! 6 EDKLC8 ,8) , EDKR C 8 * 8 ) , EBKL ( 8 , 8 ) ,EBKRC8 ,8; 
COMMON EMOD,EINER,EAREAtELEM,ERHO 
common EELK, EDKL, EDKR, EBKL, EBKR 
COMMON DW,DRAD,DLEN 

COOMON SXXL,SXYL,SYXL,SYYL,DXXL,DXYr ,DYXL,DYYL 

common sxxr,sxyr,syxp,syyr,dxxr,dxyr,dyxr,dyyr 

DMASS=DWCJ+1 )/386.4D0 

TINER=CDMASS/12,0D0)*(3.0D0*DRADCJ+n**2+DLENCJ+l; 
PINER = CDM4SS/2,0D0) *CDRADCJ + n**2) 

CALL ZER0MCCEDKR,8 .8) 

EDKRC5, 5)sDCMPLXC-DMASS*R0MEGA**2,0.0) 
EDKPC7,7)=EDKRCS,5) 

EDKRC6,6)=OCMPLX(-TINER*RDMEGA**2,0.0) 

EDKRC8,S!=EOKRC6,6) 

FDKRC6 ,8)=DCMPLX(0.0 ,-PINER*R0MEGA**2) 

EDKRC8 ,6)=DCMPLXCO.O,PINEP*ROMEGA**2) 

RETURN 

END 


OFTGINAL FAi;f "F 
OF wr, . 




,DCMPLX 


,DCMPLX 


) 




C LEFT P-NP WEARING EFFECTS RXP CCNPLEX DOUFLE PRECESION 

SUPRnUTINE EBEARL(ROMEGA) 

IMPLICIT PEAL>t=B(A-H,D-ZJ 
REAL=5‘8 EINB:RC7) .EAREAC7) ,ELFMC7) 

REAL*8 EELKCB,8) .DWCfl) ,DRAD(8) ,DLEMf8) ^ ^ 

COMPLEXi«l6 EDKL(8,8) ,EDKR(8,8) , EBKL ( 8 , 8 ) , EBKR ( 8 , 8 ) ,DCMPLX 
COMMON EMOD,ETNER,EAREA,ELEN ,ERHO 
common EELK.EDKL.EDKR, EBKL, EBKR 
common DW,DRA0,DLEN 

common sxxl,sxyl,syxl,syyl,dxxl,dxyl,dyxl,dyyl 
common sxxr,sxyr,syxr,syyr,dxxr,d::yr,dyxr,dyyr 

CALL ZEP0MCCEBKL,8 ,8) 

EBKLCl , l)=DCMPLXCSXXL,ROMEGA*OXXL) 

EBKL ( 1 , 3)=DCMPLX(SXYL,RUMEGA*DXYL) 


FBKLCi 

EbKL(3 

FRKLC3 

RBJTUPN 

END 


n=DCMPLX(wSYXL,ROMEGA*PYXLJ 

3)=nCMPLXCSYYL,ROMEGA*DYYL) 


C RIGHT END BEARING EFFECTS 8X8 COMPLEX DOUBLE PRECESIQM 

SUBROUTINE EBEA RR C ROMEGA ) 

IMPLICIT REAL*8 CA-H,0-ZD 
C RIGHT END BEARING EFFECTS 8X8 COMPLEX 

REALMS EIMERC7) .EAREAC7) .ELENC7) 

REALMS EELKC8 ,8) ,DWC8) ,DRADC8) ,DLEN(8) 

complex* 16 EDKLC8,8) , EDKR ( 8 , 8 1 , EBKL ( 8 , 8 ) , EBKR C 8, 8) ,DCMPLX 
common aMOD,EINER,EAREA,ELEN,ERHO 
common EELK,EDKL,EDKR, EBKL, EBKR 
common DW,PRAD,DLEN 

common SXXL,SXYL,SYXL,SYYL,DXXL,DXYL,DYXL,DYYL 
common SXXR,SXYR,SYXR,SYYR,DXXR,DXYR,DYXR,DYYR 
CAT'L ZSR0MCCEBKR,8,8) 

E8KR(5,5)=DCMPLXCSXXR ,R0MEGA*DXXR) 

E6KRC5.7 =DCMPLXCSXYR,ROMEGA*DXYR) 

EBKRC7 ,5)=0CMPLX(SYXR,R0MEGA*DYXR) 

EBKRC7 ,7)=DCMPLX(SYYR,R0MEGA*DYYR) 

RETURN 

END 





C LPTT Ff-l' ni«^K tiFFECTS 8Xt^ COf'Pt.FX COUHLR n^RCli^TON 

-S ii f > W' n n rj t K n t s K L ( J , R 0 M f G A ) 

TV'PLTCIT fRAL«8 CA-H,0-Z) 

, »EAr,,«y FINEHC?) ,EAREA(7) ,ELF:N(7) 

RRAL*8 EFI.»^C8,8) .DWC8) ,DRAnC8) ,DLEN(H) 

COMPt.,£X#U EDKL(8,R) ,EDKPC8,R) ,EBKI,(8,8) ,b;8KRC8,8 

cqmmhn emod,einer,earea,elen .ERHO 
COMMON- eeuk,edkl,eokr,ebkl,ebkr 

’ COMMn\- DW,DRA[),DLEN 

common SXXL,SXYL,SYXL,vSYYLiDXXL,DXyL,DyXL,DYYL 
CqMMO.N^SXXR,SXYR,SYXB,SYYR,DXXR,DXYR,DYXR,DyYR 
nMASS=DW(J3/386.4DO 

TINEP = CDMASS/l2.0DO3»!‘C3.ODO!«DRADCJ)^‘'*i'2 + DLEN(J3>!'*2 
PINEP=CDMAS5/2.0D0)*(DRADCJ3^^2) 

CAr.L ZERQMCCEDKL,8 ,8) 

fDFJjC I , nsDCMPLX(-DMASS=f<ROMEGA#-1^2,0.0) 

^OFLC3. 33=EDKLC I , n 

EDKL('2,2) = OCMPLXC“TIMER*ROMEGA4:>X2,0.0) 
EDFL(4,4)=EDKL(2,23 

FDFL ^ 2 , 4 ) sDCMPLX C 0 , 0 , -PINER>«‘R0MEGA*>I<2 3 
FnKLC4,23=OCMPOXCO.O ,PIMER^ROMEGA*^23 
RETURN 
END 


8: 


C RIGHT END DISK EFFECTS 8X8 COMPLEX DOUBLE PRECESION 

SUBROUTINE ED ISKR ( J , RQMEGA 3 
implicit REAL=t<8 CA-H.O-Z3 
REAL*8 EIMEPC73 .EAREAC73 .ELENC73 
PEAl,*8 EELK(8,83 ,DWC8 3 .DRADC8 3 ,DLEN(83 
COMPLEX’!' 16 EDKLC8,8 3 , EDKR C 8 , 8 3 , EBKL C 8 , 8 3 ,EBKRC8, 
common EMQD^EINER.EAREA ,elem ,ERH0 
EELK ,EDKL,E:DKP/EBKL,EBKR 
CW,DRAD,DLEN 

SXXL,SXYL,SYXL,SYYL.DXXL,DXYL,DYXL,DYYL 
SXXR,SXYR,SYXR,SYYR,DXXR,DXYR,DYXR,DYYR 
nMASS=DW(J+l 3/3R6.4D0 

timep=cdmass/i2.odo3=»‘C3.odo’!'Dradcj+i3*«2+dlencj+i ; 
PlNER = CDMftSS/2 .0D03’!‘CDRADCJ + 1 )’!‘*23 
CALL ZEROMCCEDKR,8 .83 
EDKRC5, 5 3=DCMPLXC-DMASS*ROMEGA’!'=!'2,0.0 3 
EDKPC7,73=EDKRC5,53 

RDKRC6,6 3=OCMPLX(-TINER*ROMEGA’I'>!'2,0,0 3 
EDKRr8,S3=EOKR{'6.63 

FOKR C 6 , 8 3 =DCMPLX C 0 . 0 , "PINER*ROMEGA’!'’!‘2 3 
EDFRC8 ,63=DCMPLXCO.O ,PINER*R0MEGA=«>!'23 
RETURN 
END 


common 

COMMON 

COMMON 

common 




3 ,DCMPLX 


,DCMPLX 


2 3 




C LEFT END REARING EFFECTwS 8X« CCMPLEX DOUBLE PRECESION 
SUBROUTINE EBEARL ( ROMEGA ) 


IMPLICIT PEAL»8(A-H,0-Z) 

REAL>!‘8 EINERC7) .EAREAC7) ,ELENC7) 


nc.nM'^P 

REAL*8 EELK(8»8) tDW(fi) ,DRAD(8) ,DLEMf8) 

COMPLEX=»‘16 EDKLC^,8) ,EDKRC8,8) ,EBKL(8,8) ,EBKR(8,8) ,DCMPLX 
COMMON EMOD,ETNER,EAREA,ELEN,ERHO 
common EELK,EDKL,EDKR,EBKL,EPKR 


common DW,DRAD,DLEN 

COMMON SXXL,SXYL,SYXL,SYYL,DXXL,DXYL,DYXL,DYYL 
COMMON SXXR tSXYR,SYXRfSYYR,DXXR,DXYP,DYXR ,DYYR 
CALL ZEPOMC(EBKL,8,8) 

EBKLCl , n=DCMPLXCSXXL,ROMEGA*OXXL) 

EBKLC i , 3)=DCMPLXCSXYL,R0MEGA»DXYL) 

EbKL(3, n=DCMPLXCSYXL,RQMEGA*DYXL) 

FBKLO , 3)=nCMPLX(SYYL,ROMEGA*DYYL) 


RETURN 

END 


C RIGHT END BEARING EFFECTS 8X8 COMPLEX DOUBLE PRECESIOM 

SUBROUTINE EBEARR (ROMEGA ) 

IMPLICIT REAL=«‘8 (A-H,0~Z) 

C RIGHT END BEARING EFFECTS 8X8 COMPLEX 

REAL>*«8 EIMERC7) .EAREAC7) .ELENC7) 

REALMS EELKC8,8) ,DWC8) , DR AD (8) ,DLEN(8) 

complex* 16 EDKLC8,8) , EDKR C 8 - 8 ) , EBKL ( 8 , 8 ) ,EBKR(8,8) ,DCMPLX 
common aMOD,EINER,EAREA,ELEN,ERHO 
common EELK,EDKL,EDKR,EBKL,EBKR 
COMMON DW,nRAD,DLEM 

common vSXXL,SXYL,SYXL,SYYL,DXXL,DXYL,DYXL,DYYL 
common SXXR ,SXYR,SYXR,SYYR,DXXR,DXYR,DYXR,DYYR 
CAt'L ZEROMCCEBKR,8,8) 

EBKR(5,5)=DCMPLXCSXXR ,R0MEGA*DXXR) 


DYYR 


EBKRC5,7)=DCMPLXCSXYR» R0MEGA*DXYR3 
EBKRC7 ,51=0CMPLXCSYXR,R0MEGA*DYXR) 


EBKRC7 ,71=DCMPLXCSYYR,R0MEGA*0YYR) 

RETURN 

END 






.SURRPUTIME ZEROMCA.IrJ) 

C X^ITiaLTZES A REAL MATRIX TO ZE^Q 
PF,AL=«=8 Ad) 

TI=I*J 

nc 10 K=l,It 
10 aCK)=0.0D0 
RETURN 
END 


C 


SUBROUTINE ZEROMC ( A , I , J ) 
INITIALIZES A COMPLEX MATRIX TO ZERO 
COMPLEX’Kie ACI.J) 

DO 10 K=1,I 
no 10 L = 1,.J 

10 ACK,r,) = (O.ODO,O.ODO) 

RETURN 

ENn 


SUBROUTINE CADDMl ( A , B , I , J ) 

C ADDS A REAL MATRIX TO A COMPLEX MATRIX 
REAL>i<8 A(I,J) 

C0MPLEX*16 BCI,J) 

DO 10 K=1,I 
DO 10 L=1,J 

10 BCK,L)=ACK.L)+SCK,L) 

100 RETURN 
END 


SUBROUTINE CMULTM C A , K , I , 8 , J , C ) 
multiplies two COMPLEX MATRICES 

COMpLEpio A (32,32 ),BC 32,1), C( 32,1) 

DQ 10 

DO 10 M=1,J 
CCL,M)=(0.0D0,0,0D0) 

DO 10 W=1,I 

C(L,M)=C(L,M)+ACL,M)*BCN,M) 

10 continue 

100 RETURN 
END 


subroutine CADDM2CA,8,T, J) 
C ADDS TWO complex MATRICES 

COMPLEX^ie ACI,J) ,B(I,J) 

DO 10 K=1,I 
DO 10 L=1,J 

10 0CK,L)=ACK,L)+R(K,L) 

100 RETURN 
END 





non noo non ono 


SUPRmjTIME Cr ^: V (^^, A , AX ) 

C CT\'V. . .FINDS THE INVERSE np THE COMPLEX DOURLE 
C PPECESXnw matrix S RY OANSS ELIMINATION WITH 
C PARTI 4 L PI'^ntlMG. TWE INVERSE OF MATRIX S IS 
C STORf.n t ; mI 

C ^.s npnpc np ’.'ATRIX TO RP' TNVPJRTFD 

Tv.priCT'^ C'^WPLFX'^16 CA-N,n-Z) 
cpALsXh CnAWS,SMAX 

rr.Mpr F'^*lh snP,,32),AICJ2,3?.),A(32,3?) 

C 


no 10 1 = 1 , N 

no 10 J=1 ,N 
10 5(I,J)=ACI,vJ) 

INITIALIZE AI 


NM=N-1 

DO 100 1 = 1, NM 
AICI,I3=(1.0D0,0.0D0) 

DO 100 J=r.NM 

AICT,J + 1) = CO.,ODO,O.ODO) 

AlCj + 1 ,n = CO,ODO,O.ODO) 

100 continue 

A1CN,M)=C1.0D0,0.0D03 
DO 4i0 h=2,M 

SEARCH POR largest ENTRY IN (K-l)TH COLUMN OF S 


200 

210 


KMsF-1 

IMAX=KM 

SMAX=CDABSCSCKM,KM3) 
no ?J0 J=K,N 
IFCSmaX~CDABSCS(J,KM3 ) 3 
TMAX=J 

SMAX= CDABSCSCJ,KM3 3 
roMTlNUF 


200 , 210,210 


IF(IMAX-KM) 300,400,300 
SWITCH (K-13TH AND IMAXTH EQUATIONS 


300 no 310 J=KM,N 
TFJMp=SCKM,J3 
SCKM,J)=SCIMAX,U3 
310 SCIMAX ,U3=TEMP 
DO 320 J=1,N 
TEMP=A1(KM, J) 

AICKM, J3=AICIMAX,J3 
AICIMAX,J3=TEMP 
320 CONTINUE 


eliminate xck-13 from kth thru mth equations 

400 DO 420 I=K,N 

RS=S(I,KM)/SCKM,KM3 
DO 410 J=K,M 

410 SCI,J)=SCI,J3-RS=P'S(KM, J3 
DO 420 J=3 ,N 

420 AICI,J3=AICI,J)-RS*AICKM,J3 
A30 CONTINUE 


BACK SUPSTITUTE 
DO 500 1=1, N 

500 flICN, r)=AI('‘M)/S(N,N) 

DC 520 K=2,U 
MK = ^■ + 1 “K 
DO 520 J=1,’V 
nn 510 L=?,F 

510 AICNK,J3=AICNK,J)-SCNK,N + 2-L3=S'AI(N + 2-L,J3 
AICNK, J3=AI(NK,vI3/SCNK,NK) 

520 CONTINUE 
RETURN 
END 





Figure 1. System Used in Rotor - Bearing Dynamic Analysis Program 
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Figure 5. Midspan Unbalance Response of Uniform Rotor in End Bearings. 
Case B3. Plain Cylindrical Bearings, e = 0.50. 








HIOSPAN DEFLECTION RATIO 





Figure 6. Miclspan Unbalance Response of Uniform Rotor in End Bearings 
Case B5. Plain Cylindrical Bearings, e = O.^o. 
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Figure 7. Journal- Unbalance Response of Uniform Rotor in End Bearings, 
Case B7. Plain Cylindrical' Bearings, e = 0.20. 
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TOTAL NUMBER OF NODES 

3 
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2 

UNBALANCE POSITION, NODE NO. 

2. . 

BEARING ECCENTRICITY, e 

0.70 

T‘ 1 1 11 i ! 

1 1 


SPEED RATIO, 




Journal Unbalance Response of Uniform Rotor in End Bearings 
Case Bll, Plain Cylindrical Bearings, e = 0.70. 
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Figure 13. Journal Unbalance Response of Uniform Rotor in End Bearings 
Case B19. Axial Groove Bearings, e = 0.20. 
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Journal Unbalance Response of Uniform Rotor in End Bearings 
Case B23. Axial Groove Bearings, e = 0.70. 
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Figure 16. Midspan Unbalance Response of Unifonn Rotor in End Bearings 
Case B25. Partial Arc Bearings, e = 0.20. 
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Figure 17. Midspan Unbalance Response of Uniform Rotor in End Bearings. 
Case B27. Partial Arc Bearings, e = 0.50. 
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Figure 18. 
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Midspan Unbalance Response of Uniform Rotw. in,End Bearings 
Case B29. Partial Arc Bearings, e = 0^70. 
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Figure 19. Journal Unbalance Response or Uniform Rotor in End Bearings 
Case B31. Partial Arc Bearings, e = 0.20. 
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Figure 20. Journal Unbalance Response of Uniform Rotor in End Bearings. 
Case B33. Partial Arc Bearings, e = 0.50. 
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Journal Unbalance Response of Uniform Rotor in End Bear 
Case B35. Partial Arc Bearings, e = 0.70. 
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Figure 22. Midspan Unbalance Res 
Case B37. Tilting Pa 
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Figure 25. Journal Unbalance Response of Uniform Rotor in End Bearings 
Case B43. Tilting Pad Bearings, e = 0.20. 
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Figure 26. Journal Unbalance Response of Uniform Rotor in End Be,?iring; 
Case B45. Tilting Pad Bearings, e => 0.50, 
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Figure 28. Midspan Unbalance Response of Uniform Rotor in End Bearings 
Case B49. Plain Cylindrical Bearings, e = 0.30. 
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Figure 29. Quarter Span Unbalance Response of Uniform Rotor in End 
Bearings. Case B50. Plain Cylindrical Bearings; e = 










(d) Equivalent Rotor in 

Equivalent End Bearings 


Figure 30. Procedure for Finding an Equivalent Rotor 
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Figure 34. Response of Stepped Rotor. Case 2. Rigid Bearings 
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Figure 36. Equivalent Rotor Response. Case 2. Flexible Bearings 
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Figure 37. Response of Stepped Rotor. Case 3. Rigid Bearings. 
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Figure 39. Response of Equivalent Rotor. Case 3. Flexible Bearings 
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Figure 41 


. Response of Stepped Rotor. Case 4. Flexible Bearings 
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Figure 42. Response of Equivalent Rotor. Case 4. Flexiole Bearings 
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Figure 43. Response of Stepped Rotor. Case 5. Rigid Bearings. 
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Case 5. Flexible Bearings 








WLITUDE CKAJQR ELLIPSE HABIUS/ECCENIRICITI/ 


: / ' 


TOTAL NUMBER OF NODES 3 

AMPLITUDE AT MIDSPAN, NODE NO. 2 

UNBALANCE POSITION, NODE NO. 2 

BEARING ECCENTRICITY, e 0.0 


SPEED RATIO, ai/w. 


Figure 45. Response of Equivalent Rotor. Case 5. Flexible Bearings 
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BALANCE L_ 12 5 12 5 — ^ 

PLANE ' ' ~ ' PLANE 


ROTOR WEIGHT: 88 LB. BEARINGS: 4 SHOE TILTING PAD 

END DISKS; 18 LB. EACH L/D RATIO: 1.0 . 

CENTER DISK: 36 LB. CLEARANCE RATIO: 3 x lO'^IN./IN. 

DRIVE: ELECTRIC MOTOR AT END (MACHINED CLEARANCE) 

SPEED'RANGE: 24,000 RPM PAD ARC LENGTH: 80 DEGREES 

FLUID VISCOSITY; 0.65 CS (77*F) PIVOT POSITION: 44 DEGREES 
0.51 CS (I30"F) FROM LEADING EDGE 

GEOMETRIC PRELOAD: 0.5 

Figure 46. Lund-Orcutt Rotor System Details, 


Figure 47. 



Unbalance Response Amplitudes of Lund-Orcutt Rotor in Free- 
Free Mode. From Ref. [5]. 
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IRANSOUCER LOCATIONS 


Figure 50. Lund-Tonneson Rotor System Details. 



Figure 51 . 


Unbalance Response Amplitudes of Lund-Tonneson Rotor. 
From Ref. [16]. 
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Figure 53. Kendig Gas Turbine Rotor System Details. 
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